Raman scattering in a d-wave superconductor - a one gap scenario 

A. V. Chubukov^ and M. R. Norman^ 



^Department of Physics, University of Wisconsin, Madison, WI 53706 



Materials Science Dtvision, Argonne National Laboratory, Argonne, Illmots 60439 



(Dated: February 28, 2008) 

Recent Raman scattering data in Big and -829 geometries in the superconducting state of un- 
derdoped cuprates were interpreted as evidence for two distinct energy gaps. We argue that these 
data can be equally well explained within a one gap scenario if final state interactions are taken into 
account. In particular, we show that they can account for the differing doping dependences of the 
Raman peaks in these two geometries. 
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I. INTRODUCTION 

A key issue in the physics of cuprates is the relation 
between the pseudogap and the superconducting gap^. 
One class of theories implies that the pseudogap and the 
superconducting gap are of different origin. The pseu- 
dogap appears at T* in the antinodal region as a re- 
sult of some competing instability (with or without long 
range order), with a magnitude that increases with un- 
derdoping, scaling with T*. The superconducting gap, 
on the other hand, appears only below Tc on an arc 
of the Fermi surface around the node that was not re- 
moved by the pseudogap, and has a magnitude that de- 
creases with underdoping, scaling with Tc- This can be 
contrasted with another class of theories that assumes 
there is a single d— wave pairing gap (not necessarily of 
the form cos fc^a — cos fc^a), whose magnitude scales with 
T* both in the nodal and antinodal regions. Above Tc, 
superconducting coherence is destroyed by fluctuations, 
but gap-like features in the spectral function survive as 
long as the magnitude of the fermionic damping, ri{T), 
is smaller than the angle-dependent pairing self-energy, 
A(0) (specifically, |A(0)| > y^r]{T)^). This give rise to 
a T dependent arc around the node, where along the arc, 
|A((^)| is smaller than \/3ri{T), and the spectral function 
has a peak at the Fermi energy, as in the normal state. 

One gap scenarios are consistent with recent photoe- 
mission measurements of the temperature dependence of 
the Fermi arc above T^^. Measurements by the same 
group below found a simple d-wave gap of the form 
cos 2(f> with a magnitude and temperature dependence 
unrelated to Tc, but consistent with T* . Similar behavior 
has been inferred from Fourier transforms of recent scan- 
ning tunneling datai^, except for regions near the node 
where there is some evidence for a suppressed gap as in 
earlier photoemission datai. A one gap scenario is also 
consistent with point contact SIS tunneling data^. 

Recent Raman studies of underdoped Hg— based 



cuprates by Le Tacon et al^ were interpreted as evidence 
for the existence of two distinct energy scales in the un- 
derdoped cuprates. The Raman spectra in Big and i?2g 
geometries behave differently with doping: the peak in 
Big geometry shifts to a higher energy with underdoping 
and, to a first approximation, tracks (twice) the antin- 
odal gap, while the peak in i?2g geometry shows the same 
doping dependence as Tc and shifts to a lower energy 
with underdoping. As Big and i?2g responses predom- 
inantly come from the antinodal and nodal regions, re- 
spectively^*', this result was interpreted as evidence that 
the nodal and antinodal gaps have different doping de- 
pendences. This does not necessary imply two distinct 
gaps (the interpretation of Rcf. 9 assumes a rf— wave gap 
with doping dependent anisotropy as in Ref. 0), but still 
it does imply that the gap function behaves differently in 
these two momentum regions. 

Here, we argue that the Raman data can be explained 
equally well within a one gap scenario, as a result of final 
state interactions, which are assumed to increase with 
underdoping. We find that final state interactions can 
lead to different doping dependences of the Big and i?2g 
responses even if the superconducting gap has a simple 
cos 2(/) form. 



II. FORMALISM 

Without final state interactions, the Raman intensity 
due to electronic excitations is given by the Fermi surface 
average of the sum of normal (GG) and anomalous {FF) 
bubbles, weighted with Raman vertices^i. In a situation 
where the self-energy depends only on w, but not on k 
normal to the Fermi surface (apart from a trivial veloc- 
ity renormalization), the Raman intensity at T = in 
the geometry labeled by 'i' is the imaginary part of the 
Raman bubble: Ri{^l), which for > is given by 
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where i?o is a normalization factor, (I>± = lu± — E(a;±), 
w-t = a;±r2/2, (/)is the angular variable along the Fermi 
surface, < ... >fs denotes averaging over the Fermi sur- 
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for a circular Fermi 



face (< ... >_F5= (2/7r) /J 
surface), and A(0) is the pairing self-energy which for 
simplicity we assume to be independent of frequency. For 
a circular Fermi surface, we assume a d— wave gap of the 
form A(0) = A cos 20. The vertices ji are different for 
different scattering geometries and are jBig {4>) — cos 20 
(i.e., cos/cajO — cosfcj^a), and 7^23(0) — sin 20 (i.e., sinfc^ja 
sinfcyfl). Because of the angular dependences of 7i, Big 
Raman scattering predominantly probes the antinodal 
regions, where cos 20 is the largest, while the B2g Raman 
intensity comes from the nodal regions, where sin 20 is 
the largesl3iiii^. 

In the BCS approximation, E = 0, and the Raman 
intensity reduces to^^^ 
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The expressions for ImRB^ and ImRB2 in this approx- 
imation can be analytically expressed in terms of com- 
plete elliptic integrals^^. At small frequencies, only nodal 
fermions contribute to the Raman intensity, and 
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where Vn is the nodal 'velocity' (u„ — 2 A for a cos 20 
gap). The real parts of the Rb^ and i?^^ bubbles be- 
have as 
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At larger frequencies, ImRsi (f^) diverges logarithmi- 
cally at 2A, while ReRBigi^) jumps from a positive to 
a negative value: 
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The intensity {ImRB2g i^)) has a broad maximum at 
around 1.6A. The real part of the B2g bubble is negative 
for all frequencies, and has a weak minimum at 2 A, where 
ReRB,(n) = -4i?o/3. 



The imaginary and real parts of Rbi,, 
general are shown in Figs. [1] and [D 
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FIG. 1: (Color online) The imaginary part of the Raman 
bubble, B\g (solid line) and -B29 (dashed line), with f2 in 
units of 2A. In the BCS approximation, these are the B\g 
and B^g Raman intensities. The logarithmic divergence of 
the B\g intensity at 2A is smoothed by the presence of a 
small damping, r\ = 0.06A. Here and in other figures, we plot 
the Raman intensity in units of ARo/n. 
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FIG. 2: (Color online) The real part of the Raman bubble. 
Big (sofid line) and B2g (dashed line). Observe that Re-Rsig 
changes sign, but Rei?s2g remains negative for all frequencies. 

Final state interactions arise from multiple insertions 
of the fermion-fermion interaction into the Raman bub- 
ble. The fully renormalized four-fermion interaction 
ai3,-ys{<lT^) generally has components in the spin and 
charge channels, and depends on the transferred mo- 
menta q — k ~ k' and transferred frequency = uj — to' . 
Restricting to only Big and i?2g harmonics, we can ap- 
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proximate T as 



where indices c and s refer to charge and spin, and F'^^* = 
V'^icu -oj'). 

We make two assumptions in evaluating the effect of 
final state interactions. First, we assume that the random 
phase approximation (RPA) is valid, with the full Raman 
intensity given by 
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where F^ = F^ + 3F|. Second, we assume that the in- 
teractions F'^'^(r2) depend only weakly on frequency for 

<< W, where W is the fermion bandwidth, but are 
strongly reduced at frequencies comparable to W. This 
holds if, e.g., the effective interaction is mediated by over- 
damped spin fluctuations. This is relevant to final state 
interactions because the constant term in Ri (the '2' in 
the r.h.s. of Eq. ([T])) comes from fermions with energies 
comparable to W, while the rest comes from fermions 
with energies of order << W. Because the effective 
interaction is strongly reduced at fi ^ W, the constant 
term should be dropped from ReRi in the the denomina- 
tor of Eq. At the same time, for the rest of Ri, the in- 
teractions F' can be safely approximated by constants^^. 

The effect of the final state interactions obviously de- 
pends on the sign of F^. The interaction in the Big chan- 
nel is the same one as gives rise to d— wave supercon- 
ductivity. If the spin component of F dominates over 
the charge component, Fs^^ is negativ o^^'^^ . F^^^ on the 
other hand is not related to pairing, and in general can 
have either sign. 



Note that the maximum value of the gap is still A. In 
Figs. [3] and m we show the Big and i?2g Raman profiles 
for a = and a finite a versus ft. We see that the peak 
position in i?2g geometry progressively shifts to lower fre- 
quencies with increasing a, and flattens. The Big peak, 
on the other hand, doesn't move, but just broadens. Be- 
cause A increases with underdoping, the Big peak, lo- 
cated at 2 A, actually shifts to higher frequency with 
underdoping. The i?2g peak, however, shifts to lower 
frequency with underdoping if the doping dependence of 
a overshadows the growth of A. We also note that the 
slope of the B2g Raman intensity at = increases with 
increasing a. 
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FIG. 3: (Color online) Big Raman intensity for a gap 
anisotropy parameter a = (solid) and a — 0.2 (dashed), 
with 7] = 0.06A. Observe that the peak doesn't move but 
just broadens with a. 



III. RESULTS 

A. General d wave gap 

As discussed above, the peak positions in the Big and 
B2g Raman intensities as a function of doping do not 
scale with each other. One way to account for this is to 
still use the BCS formula, but assume that the gap A(0) 
progressively deviates from the simple cos 2^ form with 
underdoping, such that the gap near the node probed in 
B2g scattering, and the antinodal gap probed in Big scat- 
tering, have different doping dependences (a 'two scale' 
scenario). A simple way to account for this is to add a 
cos6(/) harmonic to the d^2_y2 gap functionSiS.: 

A((^) = A[(I-a)cos2(/)-|-acos60] (8) 



B. Final state interactions 

We now show that final state interactions also lead to 
distinct behaviors of the Big and B2g Raman intensities 
(Figs. 5-8), even for a simple d— wave gap (i.e., a = 0). 
In Fig. Owe show the Big Raman intensity for negative 
Tsig- From here on, we express F in units of 7r/(4i?o)- As 
l^Big I increases, the Big peak shifts to a somewhat lower 
frequency compared to 2A, and its intensity decreases. 
Still, the actual frequency of the Big peak will increase 
with underdoping if the increase of A overshadows the 
reduction of the peak frequency relative to 2 A. 

For F^i > 0, the effect of the final state interaction 
is opposite - the peak shifts to a higher frequency, and 
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FIG. 4: (Color online) B2g Raman intensity for a gap 
anisotropy parameter a = 0,0.1,0.15,0.2, with rj = 0.06A. 
The peak moves to lower frequency with increasing a. 



its intensity becomes quite large for F^ig ~ 1- For larger 
, the intensiy drops. We show this behavior in Fig.[6l 
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FIG. 5: (Color online) Big Raman intensity for A((^) = 
A cos 20 with final state interaction Fsig = (dashed) and 
Fflj^g = — 1 (solid). For finite Fs^^, the intensity decreases in 
magnitude, and the peak shifts to a somewhat lower frequency 
compared to 2A. Here, and in other figures, F is expressed in 
units of 7r/4i?o. 

We found a stronger effect of the final state interac- 
tion in geometry. For negative F^^gi the intensity 
decreases, and the peak frequency rapidly shifts to a 
smaller value, demonstrating the same trend as with the 
gap anisotropy parameter, a. We show this behavior in 
Fig. [71 The only qualitative difference between this figure 
and Fig. |4] is that in Fig. [71 the slope of the B2g Raman 
intensity is unaffected by the final state interaction sim- 
ply because near 17 = 0, '1' dominates over i?(ri) in the 
denominator of Eq. ([7]). 

For positive Fs^^ the effect is the opposite - the peak 
shifts to a larger frequency, and its intensity increases 
quite dramatically for F^a ~ 1 (see Fig. [5]). For larger 
Fb2 , the intensity drops, and the profile becomes almost 
flat." 
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FIG. 6: (Color online) B\g Raman intensity with the final 
state interaction of opposite sign compared to Fig. [S] The 
peak shifts to a higher frequency and its intensity increases 
for Fsjg ~ 1. For larger Fs^^, the intensity drops. 
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FIG. 7: (Color online) B2g Raman intensity with the final 
state interaction Fbj^ = 0, —0.5, —1, —1.5. With increasing 
^B-ig , the peak progressively shifts down in and the in- 
tensity decreases. The slope at S7 = 0, however, remains 
unchanged. 



1. Theoretical reasoning 

The behavior of the B^g intensity for a d— wave gap in 
the presence of final state interactions has been studied 
befor o^^d^d'^ . The reduction of the peak frequency to be- 
low 2 A for negative F^j^ was understood as an exciton- 
like effect, similar to the one which gives rise to a peak in 
the dynamic spin susceptibility in a d,j.i_yi superconduc- 
tor. ReRsigiS^) is positive below 2A and is a monotonic 
function of fl in this range (Fig. The denominator 
of Eq. ([7]) contains 1 -|- TBigReRsigi^), and for strong 
enough, negative Tsig (from Eq. (O, {Tsigl > 3/7r), 
1 -I- TBigReRBigi^) equals zero somewhere below 2A. 
In the absence of damping, this would imply a pole, the 
same as for the spin susceptibility. For the Raman inten- 
sity, however, ImRBigi^) is non-zero for all frequencies. 
Still, ImRBigi^) scales as fl^ and is generally small be- 
low 2A. The net result is that the '2A' peak shifts to 
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FIG. 8: (Color online) B2g Raman intensity with the final 
state interaction of opposite sign compared to Fig. [T] The 
peak shifts to a higher frequency and its intensity increases 
for Vs^g ^ 1 (the same effect as in Fig. O. For larger rs2g, 
the intensity drops. 



a lower frequency with increasing jr^^ j^liiS. For large 
enough IFb^^I, the peak position and the peak intensity 

scale as y^—l/T, as can be easily derived from Eqs. ([3]), 
(|31), and ([7]). We note that the intensity of the peak 
decreases compared to the case with no final state inter- 
actions. For l^Bigl < 3/7r, a pseudo-resonance does not 
occur. Still, we find the peak position and peak intensity 
monotonically decrease with increasing iFsj^l. 

For the opposite sign of Tb^^ (Tsig > 0), there also 
exists a pseudo-resonance, but this time for fl > 2A 
(Fig. [HI). The reason is that ReRs^gi^) changes sign 
and becomes negative above 2 A (Fig. [2]). For positive 
Tsig ^ 1, 1 + ^BigReRsigi^) crossBS zero at some ^ 
above 2A, for which ImRBigi^) again small (Fig. [T]). 
For larger positive Tsig , a pseudo-resonance does not de- 
velop, and the intensity drops because of the presence of 
F^^ in the denominator of Eq. ([7]). 

For the i?2g case, ReRs^gi^) at small frequencies has 
the opposite sign (Eq. ([4])). Applying the same logic 
as above, one would then expect a shift of the B2g 
peak to a lower frequency for F^^g > 0, as in this case 
1 -|- T B^gReRB2g{^) crosses zero below 2A. However, 
Figs. [7] and [8] show the opposite trend - the peak shifts to 
a higher frequency for Fs^ '^-^ 1 , and to a lower frequency 
for T B2 < 0- The explanation is that at small frequen- 
cies, the imaginary part of RB2g (^) scales linearly with 
57, and is much larger than ImRBigi^) oc fl^. Because 
ImRB2 {^) is large, Eq. ([7]) is not enhanced even when 
^B2gR^RB2g{^) = ^1- this situation, the trend with 
positive TB2g is determined by the fact that ReRB2g (^) 
doesn't change sign with - it passes through a weak 
minimum at = 2A and remains negative at higher fre- 
quencies. For TB2g ^ 1, the B2g intensity has a pseudo- 
resonance at il > 2A, just like the Big intensity does. 
As the imaginary part of the Raman bubble ImRB2g i^) 
is small above 2 A (Fig. [T|), the intensity of this pseudo- 



resonance is large, and it dominates the profile of the 
Raman intensity for TB2g ^ 1- For larger Fs^^ , a pseudo- 
resonance does not develop ( \rB2gReRB2gi^)\ > 1)7 and 
the intensity drops, similar to the Big intensity. 

For negative Fsj^ , a pseudo- resonance does not occur, 
and the trend with Fs^ < is determined by the fact 
that the final state interaction monotonically decreases 
the intensity of the i?2g Raman response as the frequency 
increases. At the smallest ft, this effect is vanishingly 
small, but for H. ~ 2A, the intensity is reduced quite 
substantially. As a result, the B2g Raman intensity de- 
velops a maximum at a frequency which becomes smaller 
as |Fs2g| increases, and the intensity at the maximum 
progressively decreases, as shown in Fig. [71 



C. Fermionic self-energy 

For completeness, we also analyzed the effect of the 
fermionic self-energy on the Raman profile. As in previ- 
ous work^i^, we assumed that the self-energy T,{uj) can 
be approximated by S(a;) — —irisgii(uj) . In the superco- 
ducting state, rj is small, but it rapidly increases above 
Tc and can affect the Raman profile in the pseudogap 
regime. 

The results for the cos 20 gap are shown in Figs. [9] 
andfTUl The peaks in both geometries get broader when 
damping increases. In addition, the Big intensity in- 
creases at small with increasing 77, while the B2g in- 
tensity decreases at small il. We verified that the same 
trend presists for a general d-wave gap (i.e., a non-zero), 
and when final state interactions are included. 
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FIG. 9: (Color online) Big Raman intensity for 77 = 0.06A 
(dashed) and rj = 0.8 A (solid). With increasing rj, the peak 
broadens, and the intensity at small frequencies increases. 



1. Universal slope for B'zg intensity 

We note in passing that a finite 77 gives rise to a new 
universal regime in the i?2g response. Namely, for small 
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FIG. 10: (Color online) B2g Raman intensity for rj — 0.06A 
(dashed) and ri — 0.8 A (solid). The peak broadens and the 
slope at n = decreases. 



< 77 << A, the B2g Raman intensity is linear in fi and 
depends only on the nodal velocity, but is different from 
that in Eq. ©. 

ImRB,^m = -Ro — , ^<V (9) 

This universality of ImRsig i^) has the same origin as 
the universal conductivity of a d— wave superconductor—. 
We found, however, that this universal behavior sets in 
only for small i] < 0.1 A, while at larger ?/, the B2g Raman 
profile at small frequencies is almost independent of the 
nodal velocity. 



IV. COMPARISON WITH EXPERIMENT 

We demonstrated that there are two possible ways to 
explain the Raman data. One explanation assumes that 
the gap anisotropy changes with underdoping. Another 
explanation is that the distinct behavior of the Big and 
B2g peaks is the effect of final state interactions (Figs. [S] 
and [7]). In both cases, one can reproduce the doping de- 
pendences of the Big and B2g peaks. The second expla- 
nation is consistent with recent photoemission data that 
finds a simple d-wave gap^, even in underdoped samples. 
The qualitative difference between the two scenarios is 
the behavior of the slope of the B2g Raman response ~ 
the slope, measured in units of Rq/A, varies with doping 
via the varying gap anisotropy, but is doping independent 
if the effect is due to final state interactions. 

To extract the doping dependence of the slope from 
the data, one needs to normalize the theoretical results 
in the same way as in experiment. Le Tacon et al.^ nor- 
malized the measured Raman response for all dopings by 
imposing an empirical 'sum rule' 

j dnn imRB^,{^) = Cp (10) 



where p is the doping and C is a constant. They argued 
that this sum rule can be approximately derived for a 
weakly doped Mott insulator. With this normalization, 
the slope of /mi^Sj is essentially doping independent, 
although the peak positions in both the Big and i?2g 
channels vary substantially. 

The theoretical Raman intensity contains a prefactor 
i?o (see Eq. (P). It is proportional to the fermionic 
density of states, A^, and the square of the quasiparti- 
cle residue, Z. Both N and Z are weakly dependent on 
doping for overdoped samples, but become doping depen- 
dent in the strong coupling regime. The normalization 
used by Le Tacon et al. implies that 

RoAl J dzzImRs^^iz) ^ Cp (11) 

where z = fi/A, and ImR{z) is the function plotted in 
the figures. 

The slope of ImRB2g{z) scales as i?o/A if the gap has 
a simple d— wave form, even if final state interactions 
are relevant (the final state interaction does not affect 
the slope). The doping independence of the slope then 
implies that Rq ^ A. Substituting this into Eq. (fTTjl . we 
find that this is consistent with the data if 

1(^32^ = J dzzImRB2g[z) (X ^ (12) 

At the doping dependence of A is known, Eq. ([72]) 
is a parameter-free condition that can confirm or dis- 
prove our theory. We computed the integral of ImRB2 
from Fig. Hand found /(O) « 0.94, /(-I) w 0.30, and 
/(— 1.5) « 0.20. Extracting the peak position in units 
of A from Fig. [7] and A from Fig. 2 of Ref. i, we 
find that F^ag = roughly describes the optimal doped 
Tc = 95K sample {p — 0.15), F^a, = — 1 the underdoped 
Tc = 86K sample {p = 0.12) for which A is about 1.4 
times larger than at optimal doping, and F^^g = — 1-5 the 
underdoped Tc = 63 K sample (p = 0.1) for which A is 
about 1.5 times larger than at optimal doping. Compared 
to optimal doping, p/A^ in the two underdoped cases 
is reduced by 0.29 and 0.20, respectively. These num- 
bers are in good agreement with our /(— 1)/J(0) = 0.32 
and /(— 1.5)//(0) = 0.21. This good agreement is a 
strong argument for an explanation based on final state 
interactions rather than for a strong doping-dependent 
anisotropy of the gap. 

It is instructive to further compare our explanation 
of the data with the one presented by Le Tacon et 
al.^ . They included both a cos(60) contribution to the 
gap A(0) (as in Eq. ([8])) and a Fermi liquid renormal- 
ization. A, of the Raman vertex. They argued that 
the vertex renormalization factor A is angle dependent, 
A = Ao(l — C cos^ 20), and that C increases with under- 
doping. For nodal fermions {(j) = 7r/4), the vertex equals 
Aq and is doping independent, but for antinodal fermions, 
its values goes down with underdoping. Le Tacon et 
al. argued that they needed the C term to explain the 
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loss of intensity of the Big peak, and also to obtain a 
maximum in the response, which in the absence of 
the C term is flat. 

In our theory, the intensity of the the Big peak shifts 
down because of the final state interaction, which is 
equivalent to the renormalization of the Raman vertex 
The suppression in our case is the consequence of the 
fact that the intensity of the peak scales as y^—l/T. In 
this, we agree with Le Tacon et al. At the same time, 
the C term in the analysis by Le Tacon et al. reduces the 
Big vertex for all frequencies, while in our theory, the Big 
Raman vertex at the smallest frequencies is enhanced by 
final state interactions if ^Big < (in Fig. [SI the Big 
intensity at small Q is larger than the one without final 
state interactions). 

In another distinction from Le Tacon et al, the Big 
peak in our theory shifts down compared to 2 A. How- 
ever, the energy of the peak still increases with under- 
doping because A increase a^°'^^i^^ . 

For i?2g scattering, the effect of the C term in their con- 
sideration is consistent with our analysis. Namely, their 
vertex becomes progessively smaller with the deviation 
of a typical angle (jj from 7r/4. This obviously happens 
as the frequency increases. In our case, the ratio of the 
renormalized vertex to the bare one also progressively de- 
creases as the frequency increases, if TB2g < 0. From this 
perspective, the angular dependence of their C term also 
mimics the effect of our final state interaction in the B2g 
channel for a negative sign of ^B2g- However, in distinc- 
tion to Le Tacon et al., we argue that this effect alone 
explains the data, i.e., there is no need to invoke a doping 
dependent change of the anisotropy of the d— wave gap. 

To be more quantitative, LaTacon et al. used the ex- 
perimental A and three other doping dependent parame- 
ters: ZAq, C, and the gap anisotropy parameter, a. Two 
out of three of these parameters are fixed by (i) the 'sum 
rule' requirement, and (ii) the experimental fact that the 
slope of the normalized Raman B2g response is doping 
independent. The one free doping dependent parameter 
is chosen to fit the position of the B2g peak. 

In our theory, we assume a simple d— wave gap and cal- 
culated the vertex renormalization based on the RPA^^. 
We therefore have only two parameters: Rq, and TB2g to 
fit the same three sets of data. As we described above, we 
found that the agreement with the data is nearly perfect 
in the sense that once we fix Rq and to match the 
slope and the peak position in the B2g channel, we find 
that the experimental normalization condition is satis- 
fied, despite the strong doping variation of A. 

We also note that the dependence on the fermionic 
damping in Figs.lQlandfTOlis consistent with the data if we 
assume, like in earlier work^i^, that rj increases strongly 
above Tc. Namely, the maximum in the Big intensity 
measured by Le Tacon et al. not only shifts, but also 
broadens with underdoping, and at small f2 the intensity 
above Tc overshoots the intensity in the superconducting 
state. For B2g scattering, the intensity in the supercon- 
ducting state overshoots the intensity in the normal state 



in a broad range of frequencies, particularly for the most 
underdoped Tc = 63K sample (see Fig. 1 of Ref. 0). This 
behavior is in agreement with Fig. [TUl 

An issue left in our analysis and in the analysis by 
LeTacon et al. is the justification of the sign of the effec- 
tive interaction in the B2g channel (this is the sign of Fs^ 
in our analysis, and the sign of the C term in the analysis 
of LeTacon et al.). To fit the data, we need F^ag < 0. 
As we said in Sec. II, the full F^j has charge and spin 
components: Fs^^ = ^B2g ~^'^^Bog- spin component 
F^^^ has opposite sign compared to F^^_^ as the interac- 
tion is peaked at or near the antiferromagnetic momen- 
tum Q = (tt, tt). That is, F| cx ji{k)ji{p)xs{k—p), where 
k and p are on the Fermi surface, and the sign difference 
between the B2g and Big channels is due to the fact that 
iBigik + Q) = - lBig{k), while ■^B2g{k + Q) = lB2gik). 
If d— wave pairing is magnetically mediated, F^^^ < 0, 

hence F^^^ > (from a pairing perspective, the lat- 
ter is repulsive). However, the dominant contribution to 
T s conies from regions near the hot spots, which in the 
cuprates are rather close to (tt, 0), for which "fB2g = 0- 
As a consequence, the magnitude of F^^ is substantially 
smaller than that of iTs^gl- In this situation, the charge 
component of the i?2g Raman vertex may well exceed the 
spin component. The charge and spin components of F^ 
contribute with a different sign to the pairing channel 
(because of the spin factor cr" in the pairing vertex). Ac- 
cordingly, the charge-dominated repulsive pairing inter- 
action in the i?2g channel corresponds to < 0, and, 
hence, to a negative Fb^^, which we need. In this con- 
text, the increase in iF^^g | with decreasing p well may be 
a consequence of the reduced screening when approach- 
ing the Mott transition. 



V. CONCLUSIONS 

In summary, we argued in this paper that the data 
on Big and B2g Raman scattering by Le Tacon et al^ 
can be explained as an effect of final state interactions, 
without the need to invoke a doping-dependent change 
in the anisotropy of the d-wave pairing gap. This work 
demonstrates that care must be taken when comparing 
energy gaps derived from two particle spectroscopies like 
Raman from those obtained directly from 'single particle' 
probes such as photoemission. 
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